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Within the framework of the instanton approach we present analytical results for the following 
model problems: (i) particle penetration through a parabolic potential barrier, where the instanton 
solution practically coincides with the exact (quantum) one; (ii) descriptions of highly excited states 
in two types of anharmonic potentials: double- well X 4 , and decay X 3 . For the former case the 
instanton method reproduces accurately not only single well and double-well quantization but as 
well a crossover region (in the contrast with the standard WKB approach which fails to describe 
the crossover behavior), and for the latter case the instanton method allows to study resonance 
broadening and collapse phenomena. We investigate also resonance tunneling, playing a relevant 
role in many semiconducting devices. We show that in a broad region of energies the instanton 
approach gives exact (quantum) results. Applications of the method and of the results may concern 
the various systems in physics, chemistry and biology exhibiting double level behavior and resonance 
tunneling. 

PACS numbers: PACS numbers: 05.45.-a, 05.45.Gg. 



I. INTRODUCTION 

Semiclassical mechanics has a long history. Surprisingly however some longstanding problems still exist in the theory. 
One of them, how to describe with a sufficient accuracy states near a potential barrier top, is the subject of our paper. 
It is known that the commonly used WKB method (phase integral approach) [1] , [2] is reduced to matching of wave 
functions for classically allowed and forbidden regions. The procedure technically works for linear (or first order) 
turning points, and can be relatively simply performed only for one dimensional problems. However Id problems 
are not of great physical importance, not only since reduced dimensionality which does not allow to model many 
relevant experimental situations but as well, at least partially, since Id quantum mechanical problems can be rather 
easily solved numerically. Unfortunately efficency and accuracy of direct quantum mechanical numerical methods is 
rapidly degraded for multidimensional systems, possessing many degrees of freedom, because of the extraordinary 
amounts of computational work required to perform these studies. Furthermore, an extention of the WKB procedure 
to multidimensional systems encounters fundamental difficulties because of still unsolved matching problem for the 
multidimensional WKB solutions which become singular on caustic lines separating manifolds in phase space with 
real and imaginary momenta for each among N coordinates. Since the amount of these domains increases like Nl, 
it is a tremendous task for N > 2. And after many decades of efforts, a complete and unifying descriptions of 
multidimensional WKB solutions is still not available. 

The problem was first addresses a long ago, and some attempts to overcome the difficulties of the WKB approach 
and to improve the accuracy of the method have been performed quiet successfully. Note for example [3] where the 
authors have included into the standard WKB method additionally a special type of trajectories on the complex phase 
plane, along which the semiclassical motion is described by the Weber functions (see also [4]). However the choice of 
these additional special trajectories (which one has to include to improve the accuracy of the WKB method near the 
barrier top) depends on the detail form of the potential far from the top, and therefore for each particular case the 
non-universal procedure should be perform from the very beginning (see also more recent publications [5] , where the 
authors use some distortion of Stokes diagrams, or [6], where time dependent quantum mechanical calculations for 
anharmonic and double- well oscillators have been performed). 

Thus evidently there is some need for different from WKB semiclassical approach. One of the alternative to WKB 
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semiclassical formalism is so-called extreme tunneling trajectory or instanton [7], [8], [9], could be very effective 
for calculating globally uniform (i.e. without any singularity) wave function of the ground state. It allows to find 
semiclassical wave functions for a very wide class of potentials with arbitrary combinations of the first and of the 
second order turning points. The method recently was adapted as well for the description of low - energy excited states 
[10, 11]. One of the main advantage of the instanton approach that it can be readily extended for multidimensional 
systems using perturbative techniques (see [12] and references therein). 

However before to investigate multidimensional problems, one has to study Id potentials and Id problems which 
cannot be solved accurately by the standard WKB method. It is our main concern in the paper. The generalization 
of the instanton procedure for highly excited states is not straightforward at all, and required additional analyses. 
We consider only few relatively simple examples, but such analyses are useful for gaining insight into more complex 
systems for which even approximate theoretical results are not available. 

In fact for many interesting physical problems high accuracy calculations are out of reach by the standard WKB 
method but as we will see in a little while the instanton approach is afford to overcome the difficulties of the WKB 
procedure. Since this fact has largely gone unnoticed in the previous studies, we found it worthwhile to present the 
investigation of a few simple examples in a short and explicit form, and also to point out its practical usability. And 
besides the study (apart from the aim to illustrate the efficency of the instanton approach) is a prerequisite for an 
explanation and successful description of many relevant physical phenomena (e.g. low-temperature quantum kinetics 
of phase transitions, see e.g. [13]) where an active (reaction) path is confined effectively to one dimension. 

All examples considered in our paper, related to the fundamental problems of chemical dynamics and molecular 
spectroscopy (see e.g. the monography [9] and references herein). Symmetric or slightly asymmetric double- well 
potentials are characteristic for molecules and Van der Waals complexes with more than one stable configurations 
[14], [15], [16], [17]. The states of such systems close to the barrier top (theoretically described by the instanton 
approach in our paper) are most relevant for radiationless evolution of highly excited states. These states have a 
double nature (localized - delocalized) and the nature manifests itself in the form of wave functions which contain 
simultaneously the both components: localized in one from the wells, and delocalized between the both of the wells. 
The states close to the barrier top of decay potentials govern of thermally activated over-barrier transition amplitudes. 
For the low energy states the main reduction factor is the tunneling exponent, while the contribution of the highly 
excited states is limited by the Boltzmann factor. Our instanton calculations demonstrate that there is no sharp 
boundary between quasistationary and delocalized states. Recently two of us (V.B. and E.K.) [18] investigated the 
eigenstates of the highly asymmetric double well potential. We have shown that quantum irreversibility phenomena 
occur when the spacing between neighbouring levels of a deeper well becomes smaller than the typical transition 
matrix element. Obviously this criterion can be also applied to the states near the barrier top. Note that for the 
low energy states the asymmetry providing irreversible behavior should be very large, whereas for the states near the 
barrier top, the condition to have the ergodic behavior is not so severe, it is sufficient that the asymmetry of the 
potential is comparable to the barrier height. 

Our paper has the following structure. Section II contains basic equations of the instanton method necessary for our 
investigation. As an illustration of the method we consider a touchstone quantum mechanical problem - penetration 
of a particle through parabolic potential barrier. In this case the instanton solutions which are the asymptotics 
of the Weber equation are exact. Section III is devoted to the investigation of highly excited states in a double - 
well potential. For sake of concreteness and simplicity we study a quartic anharmonic X 4 potential. The instanton 
approach allows us to reproduce accurately not only asymptotic behavior but also a crossover region from single well 
to double wells quantization. In the section IV analogous problem for X 3 anharmonic potential is studied. Section 
V is devoted to so-called resonance tunneling phenomena, interesting not only in its own right but as well playing a 
relevant role in many semiconducting double barrier structures. Finally, in the section VI we discussed the results 
and conclude. In the Appendix to the paper we compute so-called connection matrices, which provide a very efficient 
method of finding semiclassical solutions to the Schrodinger equation in potentials having several turning points. 
It is important and significant to know the connection matrices. This is not only in itself but also for developing 
good analytical approximation. Those readers who arc not very interested in mathematical derivation can skip the 
Appendix and find all the results in the main body of the paper. 

II. PENETRATION THROUGH PARABOLIC POTENTIAL BARRIER. 

A. Instanton approach. 

Let us remind for the sake of conveniency the main ideas of the instanton approach. The first step of the approach 
derived in [7] and [8] is so-called Wick rotation of a phase space corresponding to a transformation to imaginary 
time t — ► it. At the transformation potential and kinetic energy change their signs, and Lagrangian is replaced by 
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Hamiltonian in the classical equation of motion. After this Wick rotation the standard oscillating WKB wave functions 
are transformed into decaying exponentially functions which are vanished at X — > ±00. Following to [10], [11] we will 
use slightly different formulation of the instanton method, from the very beginning assuming exponentially decaying 
real-valued wave functions. Thus taking into account that the wave functions of bound states can be chosen as real 
quantities, one can look for a solution to the Schrodinger equation in the form 



* = cxp(-7cr(X)) 



(2.1) 



where 7 is semiclassical parameter^ = mOo^o/^ where m is a mass of a particle, ao is a characteristic length of 
the problem, e.g. the tunneling distance, fl is a characteristic frequency, e.g. the oscillation frequency around the 
potential minimum, and further we will set h = 1 measuring energies in the units of frequency) which is assumed to 
be sufficiently large, and a can be called an action and for this function the first order differential equation of the 
Ricatti type should hold 







+ 7 


2dX 2 



0, 



(2.2) 



where V(X) is a potential, and e gives particle eigen states (energies). Here and below we use dimensionless variables 
(for the energy e = E/Qo, for the potential V = U/'yQo, and for the coordinate X = x/a , where E and U are 
corresponding dimensional values for the energy and for the potential). We believe that 7 > 1 and therefore cr(X) 
can be expanded into the asymptotic series 



a(X) = W{X)+ 1 - 1 W 1 {X)+ 1 - 2 W 2 {X) + 



(2.3) 



The first and the second over 7 1 order terms become identically zero, if the time independent Hamilton- Jacobi 
equation (HJE) and so-called transport equation (TE) are satisfied 



1 fdW 

2 \ dX 



V(X), 



and 



dW dA 
~dXdX 



ld 2 W 
YdJG 



A = eA. 



(2.4) 



(2.5) 



where 



A(X) = eM~Wi(X)) . (2.6) 

The essential advantage of the instanton method in comparison to a standard WKB is that in the former approach 
HJE is solved at E = 0, and therefore the classically allowed regions disappear. A price we should pay for the 
advantage is an appearance of second order turning points (unlike WKB method where there are only linear turning 
points). 

It is well known that WKB wave functions are singular at the turning points, and therefore different approximations 
represent the same wave function in various domains. The famous Stokes phcnomenom [2] is related to the distribution 
of the turning points, and Stokes and anti-Stokes lines emanate from each turning point. By the definition Stokes 
lines represent the lines where the dominance of the dominant exponential semiclassical solution to the Schrodinger 
equation becomes strongest, and anti-Stokes lines, on the other hand represent the lines across which the dominance 
and subdominance of the solutions interchange. Evidently near Stokes and anti-Stokes lines the WKB approximation 
does not work and should be refine [2]. On the contrary, since classically accessible regions do not exist within the 
instanton formalism, the Stokes lines pass continously through second order turning points, and globally uniform 
real solutions to the Schrodinger equation can be constructed using the asymptotically smooth transformation of 
the instanton wave functions into the Weber functions. Just this global uniformity is the principal advantage of the 
instanton method. 

A clearer idea of the instanton approach is obtained by the derivation of well known [1] quantization rules for the 
harmonic oscillator (V(X) = X 2 /2). For a given energy e any solution of the Schrodinger equation can be represented 
as a linear combination of the solutions of the Weber equation [22] 



(2.7) 
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where z = X^/j, and v = e — 1/2. The basic solutions of (2.7) are the parabolic cylinder functions [22], and only the 
function D v (—z) is vanished at z — > oo for argz — 0. For the argz = ir the asymptotic behavior of this function at 
z — > oo is given by [22] 

D v (-z) = z" cxp - eMin^z— 1 cxp , (2.8) 

It can vanish at z — > oo only at the poles of T(—u), and this condition gives the exact eigenvalues of the harmonic 
oscillator 

1 

e = n+-. 

Moreover since D u (—z) for integer and positive v coincide with known harmonic oscillator eigenfunctions [1], the 
instanton approach is exact for the harmonic oscillator. 



B. Tunneling through the harmonic barrier. 

For the next less trivial illustration of the instanton approach efficiency, let us apply the method to the problem of 
quantum mechanical tunneling through a parabolic potential 

mO 2 

U(x) = U -^x 2 , (2.9) 

where m is a mass of a tunneling particle, and flo is a characteristic frequency (curvature of the potential). Besides 
the potential has a characteristic space scale ao- Using Q and ao to set corresponding scales the parabolic potential 
(2.9) can be written in the following dimensionless form 

V(X) = V -^X 2 . (2.10) 

The Schrodingcr equation in these variables 

— + [ 1 2 X 2 -a 1 ]*(X)=0, (2.11) 

where 

a = 2^, (2.12) 

and introduced above semiclassical parameter 7> 1. 

The Schrodingcr equation (2.11) can be transformed into the Weber equation [22] by 7r/4 rotation in the complex 
plane X 

1 / i7T 

-zexp 



and therefore the solution of (2.11) can be represented as a linear combination of the parabolic cylinder functions D v 
[22] 

* v {z) = c l D v (z) + c 2 D v (-z), (2.13) 

where v = -(1/2) - i(a/2). 

At X — > oo only the transmitted wave exists with the amplitude (i.e. transmission coefficient) T 

v^Texp^, ( 2 - 14 ) 

while for X — > — oo one has the incident (oc exp(— ijX 2 /2)) and the reflected wave oc Rex-p(ijX 2 /2). By a standard 
quantum mechanical procedure [1] the transmission T and the reflection R coefficients can be found using known 
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asymptotics of the parabolic cylinder functions [22] at the fixed energy (i.e. at the fixed a) and it leads to a well 
known [1] expression 

; (2.15) 



1 + exp(7ra) 



Note that the solutions (2.13) are the exact solutions to the Schrodinger equation in the parabolic potential (2.10). 
Now let us apply to the same problem the instanton approach shortly described above. The solutions of HJE (2.4) 
and TE (2.5) equations which are milestones of the method can be found strait-forward and read 

X 2 

W = ±i— , A = A a X-^ 2 cxp{±ialnX/2) , (2.16) 

where the integration constant Ao determines energy dependent phases of the wave functions. From the comparison of 
(2.16) and (2.13) one can see that the instanton wave functions are the asymptotics of the parabolic cylinder functions, 
and, therefore, since the transmission T and reflection R coefficients are determined only by the asymptotic behavior, 
the values of T and R found in the frame work of the instanton approach coincide with the exact quantum mechanical 
ones at any value of the energy (of the parameter a). Let us remind that the instanton and exact quantum mechanical 
solutions for the harmonic oscillator also coincide for any energy. 

To accomplish this subsection and for the sake of the skeptical reader it is worth to mention that the WKB wave 
functions coincide with the exact solutions only at a <C — 1. In the region \a\ < 1, i.e. when characteristic size 
of the forbidden region becomes comparable with the particle wave length, specific interference phenomena between 
transmitted and reflected waves occur, and this kind of phenomena can not be reproduced in the standard WKB 
approach assuming that all turning points are independent ones. 

As an illustration in Fig. 1 we show energetical (a) dependence of the phase for the wave function reflected by the 
parabolic potential. The exact quantum mechanical and the instanton solutions (<j) in Fig. 1) are indistinguishable 
over a broad region of energies, while the WKB solution (4 l o /KB in Fig. 1) is deviated from both of them. 

C. Connection matrices. 

Our analysis can be brought into a more elegant form by introducing so-called connection matrices. In the instanton 
approach (as in any semiclassical treatment of the scattering or of the transition processes) we are only concerned with 
the asymptotic solutions and their connections on the complex coordinate plane. Thus it is important and significant 
to know the connection matrices. These connection matrices provide a very efficient method of finding semiclassical 
solutions to the Schrodinger equation in potentials having several turning points. This is a relevant starting point 
also for developing good analytical approximation. 

It is convenient to formulate the general procedure for calculating of the connection matrices for an arbitrary 
combinations of the first and of the second order turning points. After that the procedure can be applied to any 
particular problem under investigation. To do it technically, one has to extend the known for linear turning points 
procedure [2] . All necessary details of the generalization are given in the Appendix to the paper, and here we present 
only main definitions and results. For the equation 

j2 T 

— +>f q (z)*(z)=0, (2.17) 
in the semiclassical limit 7> 1 the Stokes and anti-Stokes lines are determined by the following conditions, respectively 

ReW(z) = 0, (2.18) 

and 

ImW(z) = 0, (2.19) 

where 

W(z) = f ^qjz)dz, (2.20) 

J zo 

where zq is the turning point under consideration. 
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For the harmonic potential we have only the linear turning points for real (at a > 0) and imaginary (at a < 0) 
energies. In the Appendix we calculated all the connection matrices we need. Thus to fully analyze the problem for 
the whole range of parameters, we should know only the distributions of turning points and Stokes and anti-Stokcs 
lines on the complex plane. At real turning points (a > 0) X\^ = ^(a/j) 1 ^ 2 there are 4 Stokes and 4 anti-Stokes 
lines, and two cuts on the complex plane (see Fig. 2). 

At a ^> 1 the connection matrix can be easily calculated as the direct product of the connection matrices found in 
the Appendix (M~ from (A5) and the Hermitian conjugated matrix M+) and the following diagonal shift matrix 



exp(7ra/2) 

exp(— 7ra/2) 



(2.21) 



It leads to the transmission coefficient T ~ exp(— na/2), i.e. coinciding with (2.15) in the limit a 3> 1 with the 
exponential accuracy. To improve the accuracy at smaller values of a, at the calculation of the connection matrices, 
one must take into account not only the contributions from the contours going around the turning points, but as well 
the additional contribution to the action from the closed path (with a radius 3> |-Xi,2|) going around the both points 
X\ t 2 (see Fig. 2). The procedure changes the Stokes constant T 3 (on the dashed line separating the regions 3 and 4 in 
Fig. 2), which becomes 

\T 3 \ = [l + cxp(-7ra)} 1 / 2 , 
and finally it leads to the correct transmission coefficient 



T = iT 3 _1 exp I I 



2 J ' 

which is identical to (2.15). 

In the case a < 0, the whole picture (see Fig. 3) of the Stokes and of the anti-Stokes lines and turning points, is 
turned by the angle tt/2 with respect to depicted in Fig. 2. If one bluntly takes the point X = as the low integration 
limit for the action W* (All), it leads to the following transmission coefficient 

T=l-^exp(-n\a\), 

which can be reliable (with the accuracy exp(— 27r|aQ) only for \a S> 1. Again as we did for a > case, to improve 
the accuracy we should take into account the contribution from the path surrounding the both imaginary turning 
points (this fact was noticed long ago by Pokrovskii and Khalatnikov [3]). 

At the isolated linear imaginary turning point iX\ the connection matrix is found from (A5) 

M+=(J i«p(-7r|a|/2)^ (2 22) 

Analogously the Hermitian conjugated matrix M-j~ comes from the contribution from the closed path surrounding 
—iX\. These contours provide only the amplitude of the dominant (exponentially increasing) wave. However the 
accuracy is not enough to find the amplitude of the corresponding subdominant solution (exponentially decaying 
wave function), and thus it leads wrongly to the transmission coefficient T = 1. To improve the accuracy and to find 
correctly T, one should include to the procedure the connection matrix for the isolated second order turning point 
(which is in this particular example, the maximum of the potential). Using (A9) we can find explicitely this matrix 

u _ ( (l + expHrH)) 1 ^ iexp(-7r|a|/2) \ 

-iexp(-7r|a|/2) (1 + cxp(-7r|a|) 1/2 / [ 6 > 

In principle the same kind of calculations might be performed in the adiabatic perturbation theory (which employs 
in fact the Plank constant smallness equivalent to 7 » 1). Note for example the paper [21] where the contributions 
from the contours surrounding turning points (analogous to presented above) have been taking into account. It seems 
very plausible that on this way one will be able to combine the instanton approach and the adiabatic perturbation 
theory, however, this issue is beyond the scope of our paper and will be discussed elsewhere. 



III. HIGHLY EXCITED STATES IN DOUBLE - WELL POTENTIAL 



Literally, the instanton approach described in the previous section, is valid for the states with characteristic energies 
which are small in comparison to the barrier height. However as we will show in this section the instanton method 
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works pretty good for the energy states near the barrier top V$. As an illustration let us consider symmetric double 
well potential (quartic anharmonic X* potential) 

Vo-V(X) = ±X 2 (l-X 2 ). (3.1) 

The Schrodinger equation with the potential (3.1) can be rewritten in dimcnsionlcss variables in the following form 
most convenient for the application of the instanton approach 

— + [2 7 2 (Vq - V(X)) cnM*) = o • (3-2) 
Furthermore HJE and TE instanton equations are correspondingly 

2 



1 f dW 

2 V ~dX 



= V - V(X) , (3.3) 



(3.4) 



and 

dWdA 1 U 2 W 
'dXdX + 2 \~dX 2 ~ 

Formal solutions to this set of equations (3.3 - 3.4) are even and odd instanton wave functions 

Vf=A ± (X)exp(i 1 W ± (X)), (3.5) 
where the action W± (solution of HJE) should be determined from 

^=±y/2(V -V(X)), (3.6) 
and the amplitude (pre-factor) in own turn reads as 

dW± 



1/2 / 

exp i-ia I I ^ ) dX\. (3.7) 




dX 

The quantization rules [1] are related to continuous matching of the solutions at the turning points (the second order 
turning point X = 0, and the linear turning points X = ±1 for a > and X = ±i for a < 0). The crucial 
advantage of the instanton solutions (3.5) is due to the fact that these functions have no singularities at all inside the 
barrier, since the corresponding exponents are pure imaginary ones in the classically accessible regions (unlike WKB 
solutions). Besides the general form of the instanton wave functions does not depend noticeably on whether E < Vb or 
E > Vq. This advantage allows us to include the instanton wave functions into the basis of globally uniform functions 
diagonalizing the Hamiltonian even for highly excited states. 

The described above general procedure for searching instanton solutions to the Schrodinger equation with the model 
potential (3.1) has a tricky point, what is the motivation for presenting in some details below the explicit procedure 
for the search, and new results will be emanated from our investigation. The procedure includes several steps: 

• Near the second order turning point one can use the exact solution to the Schrodinger equation (2.14) with 
ci = ±C2 for the even and odd solutions respectively. For \X\ ^> 1 from (2.14) and known asymptotics of the 
parabolic cylinder functions [22] 

exp(z/(X)) , expH/PO)" 



(3.8) 



where 

2tt 



ci 



r((3 + ia)/4) 
and 



r((l-ia)/4) r((l + ia)/4) 

exp(-^)2-W4 (27 )-i/4 ) (3.9) 



f{ X) = 1 -X*-^lnX-^ln 1 -l. (3.10) 



To have correct even and odd linear combinations conforming to (3.5) 



exp(±z/i) 

c± = Cl r ( (i ± z aA /2)/4)' (3 - n) 



where fi = (a ln-f)/4: + tt/8. 
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Near the linear turning point X = ±1, the Schrodinger equation is reduced to the Airy equation [22] 

^r-y*(y) = o, (3.12) 

where at X < 

v 2 / 3 



y = 7 

The solution vanishing at y — > oo is [22] 



X + l + ^ 

7 



(3.13) 



*(2/) = l2/|- 1/4 si n (||2/! 3/2 + J) • (3.14) 

Continuing this solution into the regions (X ± 1)V^7 3> 1 and sewing there with (3.8) we come to 

^ = exp (-i2yW + i^j , (3.15) 

where W* is the energy dependent action in the interval [X = , X = 1]. 

Comparing (3.15) and (3.11) we find the quantization rules which read for the even states as 

r((i + *)/4)_^( ...3^ (316) 

and for the odd states as 



r((l-ia)/4) V 2 

r((3 + ia)/4) 



r((3-ia)/4) 



exp(-2i 7 VF* - in) . (3.17) 



Finally from (3.16), (3.17) we come to the quantization rule which can be written in the single form for the 
both, even and odd, states 

9-vM/* 4- <MU^ - ( (W 4 ) + 2nn - arctan(t/i (na/4)) \ ( 
2jW +2<P(a)=[ /n,, UOm _ 0o • (3.18) 



(37r/4) + 2?™ - 2 arctan(i/i (7ra/4)) 

The relation (3.18) is the quantization rule we looked for, and which allows us now to pick the fruits of the instanton 
method. For the highly excited states, i.e. a <C — 1 one can find from (3.18) 

2jW* + 2<j){a) =n(n ' 



2, 

where n is an integer number. For the low-energy states a>l (3.18) reproduces the known quantization rule 

( !\ , 1 ^ ™ 
7 = n I n + o I ^ o CX P V 



where W£ is the action in the classically admissible region between linear turning point in the left well, namely 

7 Wl = jW* + cp(a) . (3.19) 

Note the essential advantage of the instanton quantization rule (3.18) with respect to the traditional WKB formalizm, 
where the quantization rules in tunneling and over-barrier regions are completely different [1]. The instanton approach 
gives the single quantization rule (3.18) which is valid in the both regions and besides fairly accurate describes the 
crossover behavior near the barrier top, where periodic orbits localized at separate wells transform into a common 
figure height orbit enclosing the both wells. 

We illustrate the results of this section in the Fig. 4, where the universal dependence of the eigen values (in the 
symmetric double - well potential (3.1) on a is plotted. For the sake of comparison we presented in the same figure 
the eigen values found by the conventional WKB procedure and by the exact quantum mechanical computation. It is 
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clear from the figure that the WKB method errors are maximal in the region of small \a\, since the oscillation period 
diverges logarithmically in this region (the particle spends infinitely long time near the second order turning point). 
On the contrary, the errors of the instanton approach are minimal near the barrier top (small \a\). 

As we already mentioned at the beginning of this section, near the potential minimum, the instanton approach 
also works very accurate. Generally speaking the instanton solutions are always correct when the deviation from the 
corresponding minimum exceeds characteristic zero - point amplitudes. Mathematically the accuracy of the instanton 
approach is based on the transformation of the semiclassical solutions into the harmonic oscillator eigenfunctions 
(what ensures as well the correct normalization of the instanton wave functions). Therefore it is naturally to expect 
that the instanton method works very accurate near the barrier top and near the potential minimum. On the contrary, 
the intermediate region, where anharmonic shape of the potential is relevant, one should expect a poor accuracy of 
the instanton method. Fortunately it turns out that the mathematical nature of the problem is on our side, and 
the instanton approach has a reasonable accuracy (of the order of the accuracy of the WKB method) even in this 
region. The fact is that the instanton wave functions are exact not only in zero but in the first order over anharmonic 
corrections to the potential approximation. It can be shown using anharmonic perturbative procedure proposed by 
Avilov and Iordanskii for WKB functions [19] and generalized in [20] for the instanton wave functions. 

Besides, what is more relevant for practical computations, the instanton wave functions (unlike WKB ones) are 
continuous near its "own" minimum. Numerical estimations show that in the intermediate energetical region, the 
accuracy of the instanton wave functions is of the order of 5 — 10%. 

To accomplish the section we present the connection matrices needed to find semiclassical solutions to the 
Schrodingcr equation in the double - well potential. Analogously to the results of the Section II, the connection 
matrix for the instanton solutions is the product of the connection matrices (A5) for the linear turning points, and 
the connection matrix for the second order turning point, which is for the case under consideration, the maximum of 
the double - well potential. Using (A9) one can find for this latter connection matrix 

/ 2(exp(7ra/2) + (1 + cxp(Tra)) 1 / 2 cos(2 7 W*)) 
\ (l + exp(7ra)) 1 / 2 sin(27VK*) 



(3.20) 



-(1 + exp(Tra)) 1 / 2 sin(2 7 Fy*) 
(l/2)(- exp(7ra/2) + (1 + exp(Tra)) 1 / 2 cos(2 7 VK*)) 

It is worth to note that the reflected wave near the barrier top acquires non trivial phase factor. The phenomenom is 
related to interference of incident, reflected and transmitted waves, and thus the phase has some geometrical meaning, 
like famous Berry phase [23]. The geometrical origin of the phase manifests itself more clearly if we remind that the 
semiclassical phase factor is determined by the probability density flow through the barrier 

t ■ T * d * 

One can look to this phase factor from a slightly different point of view, since tunneling results in the phase shift 
related to the change of eigenvalues. The quantization rules (3.18), (3.19) can be rewritten in the following form, 
which is the definition of eigenvalues e„: 

1 

£n = Tl + - + Xn , 

where n is integer numbers numerating eigenvalues, and \n is determined by an exponentially small phase shift due 
to the existence of the barrier between two wells. The phase shift Xn has the same functional form (and physical 
meaning) as the geometrical phase factor (appearing due to interference phenomena) which a quantum mechanical 
wave function acquires upon a cyclic evolution [23], [24], [25]. 



IV. DECAY POTENTIAL 



In this section we present theoretical studies of highly excited states in a decay potential, which for definitencss will 
be chosen as anharmonic X 3 potential 

V -V{X)= l -X\l-X). (4.1) 
As a first (but compulsory) step we investigate the low - lying tunneling states. 
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A. Tunneling decay of metastable states. 



We start from this simple case to pick first the low-hanging fruits, i.e. to describes states 

Vb > e n > V(X -» oo) , 



(4.2) 



what means that a local minimum is separated from continuum spectrum by a high energetical barrier, and therefore 
the quasistationary states e„ are characterized by "good" quantum numbers n. Note that a generic decay potential, 
depicted in Fig. 5, is determined by the positions of the barrier top Xo and three turning points —X\ , X = , +X 2l 
and near these points 



V(X) = { 



V a - (dV/dX) x= 



X 

2 



V - \{X-X 
(dV/dX) x=X2 (X 



Xx), \X + X!\ ^0 
-» 

'\X-X \^0 
X 2 ),\X-X 2 \^0 



> . 



(4.3) 



The potential (4.1) is just a particular example of the generic decay potential (4.3) (X\ = 1/3 , Xo = 2/3 , X 2 = 1, 
and Vb = 2/27) which we will use only for the sake of concreteness and explicit illustration, all presented below results 
are equally valid for the generic potential. As a note of caution we should also remark that in the instanton approach 
to this problem we always have deal with the only two turning points. For the low energy states the points are X 2 
and the potential minimum X = 0, for the high energy states the points are —X\ and the potential maximum Xq. 

According to (4.1) there are no turning points at X > X 2 , and at X ^> X 2 the potential can be considered as a 
constant, and therefore the wave functions asymptotically at 1 > 1 2 should coincide with plane waves. Furthermore 
near the linear turning point X = 1 the Schrodinger equation with X 3 anharmonic potential (4.1) is reduced to the 
Airy equation (3.12), whose solutions are linear combinations of the indices ±1/3 Bessel functions at real (for X < 1) 
and imaginary (for X > 1) values of the arguments 



*(w) = s/u 



B+h/3 



2u 3 / 2 



B_I_ 



1/3 



2u 3 ' 2 



(4.4) 



and 



where 



*(0 = Vc 



- 2 /*3/2 

-B+Jl/3 I 1 +B_J_ 



u = (2 7 ) 2 / 3 



/2C 3/2 \ 
1/3 ^ — 3 — J 



l-X- 



"+(1/2) 



(4.5) 



(4.6) 



at X < 1, and 



C = (2 7 ) 2/3 
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(4.7) 



at X > 1 (remind that here v = e n — 1/2). 

The coefficients B± must be chosen such a manner that at ( 3> 1 (4.5) gives the plane waves, and thus using known 
asymptotics of Bessel functions [22] 



B+ = B_cxp • 



In the classically forbidden region «» 1 the instanton solutions of HJE (2.4) and TE (2.5) matching continuously 
near the turning points the quantum mechanical solutions of the Schrodinger equation is 



*± = A± exp(±7VF) , 



(4i 



where 



a± = x- 1 ' 2 {i-xy 



1/4 



'l - yi~~x~ 

i + V^~x_ 



(4.9) 
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and 



w = ^-\(i-x)^ + A -(i-xf'\ 



(4.10) 



The wave functions of the quasistationary states, we are looking for, are the linear combinations of the instanton 
solutions (4.8) and the suitable coefficients in the combinations are determined from the condition of asymptotically 
matching to the parabolic cylinder (2.13) and Airy functions, what leads to the following equation for complex eigen 
values v 



r(- 

Since the function T(z) has a simple first order pole at z 
rate T n of the quasistationary state e„ 



(in) 

n, one can easily find the main contribution to the decay 



Q 



2 y/+l/2 2 6iH-3 



• exp 



16 , 
15' 



(4.12) 



Note that for the ground state n = (4.12) coincides to the result found by Caldeira and Legget [26]. From the other 
hand the decay rate is related to the current flow [1] at X — > +00 providing the constant amplitude of the outgoing 
wave 

-1 



En 
^0 



/ \^\ 2 dX 



d^* 

h * 

dX dX 



Introducing into (4.13) the explicit form of the wave functions (4.5) , (4.6) , (4.7) we get 

r„ 97 1 / 6 . 



O 



8 2V^ |S+|2 ' 



(4.13) 



(4.14) 



which is equivalent to (4.13). 

The decay rate according to (4.14) depends only on the normalization of the instanton wave function and on the 
amplitude of the outgoing wave. The both characteristics are determined essentially by the behavior of the instanton 
wave function in the vicinity of the turning points only. Note however that in this approximation the instanton 
computation of the decay rate (4.13) or (4.14) is satisfactory one only for the ground state, because corrections of the 
order of 7 _1 strongly increase with the quantum number n. This sins of omission can be easily relaxed if we take into 
account the X 3 anharmonic contribution into the potential as a perturbation 



r„ 
n 



2 _yiH-l/226i>+3 
7T 



(Vo 
exp (-7 



1 



576 7 



(164n 3 + 246n 2 + 1216n + 567) 



(4.15) 



The decay rate calculated according to (4.15) provides the same level of the accuracy as the WKB and exact quan- 
tum mechanical computations for 7 > 5. Out of the regime of interest, the instanton theory loses all pretence of 
predictability. 



B. Highly excited states for anharmonic X [i potential. 

In the Subsection IV A we calculated the decay rate of low-energy metastable states. For this case (the states e n 
can be characterized by the good quantum number n) period of oscillations in the well is smaller than inverse decay 
rate (e„ » r„/f2 ), and T n is determined by the probability current density flowing from the well into the classically 
admissible region (X > X for a given energy e„, see Fig. 5 and notations herein) at the condition of vanishing back- 
flow from this region to the barrier. Evidently the method does not work for highly excited states, when T n > e„ilo- 
In this Subsection we go one step further with respect to the Subsection IV A extending the instanton approach to 
the decay of highly excited states. 

First it is worth noting that the wave functions must be vanished at X — > —00, and besides could be always chosen 
as real - valued quantities at X — > +00. From these two conditions one can find the relations between the instanton 
wave functions in the regions X < X\ and X > Xq (see notations in Fig. 5), and, as a result of these relations, to 
calculate the phase 5(a) (counted from the barrier top) of the standing wave in the region X > X - It reads 

■ / -o urn l + v / 2^cxp(-™/4)exp(i27VF*)[r(l-m/2)]- 1 

exp(i2S) = -iexr>(-i2^W) = — — — — — ■ — — . (4.16) 

1 + V2^cxp(-7ra/4) cxp(-i2 7 Fy*)[r(l + ia/2)]- 1 
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According to the standard quantum mechanics [1] the phase (4.16) determines the scattering amplitude. Thus from 
(4.16) we can find the scattering amplitude in the deep classically forbidden region, and as a result of it to compute 
the eigen values in this region. To perform the calculation we need to know the terms of the order of exp(— 7r|a|) in 
the T-functions expansion (which are beyond the standard Stirling formula) [22] 

~ v^exp (-^ ± i0) (l - 1 exp(-7r|a|)) , (4.17) 



where 

" \a\ 



In- 

2 



(4.18) 



Finally from (4.16), and taking into account (4.17), (4.18) we find (with required exponential accuracy) the poles of 
the scattering amplitude 



2 7 1U£ = 2 7 IU* + <p(a) = 7r(2n + 1) - i 



7T 1 

-(|a|-a) + -exp(-7r|a| 



(4.19) 



and therefore explicitcly solving (4.19), the complex eigen values, and in particularly the decay rate for highly excited 
states in the anharmonic decay potential. 

The same manner as for low energy tunneling states, for highly excited states (i.e. at \a\ 3> 1) the real part of the 
eigen values e„ is determined by the action along closed trajectories in the well, whereas the imaginary part (i.e. the 
decay rate T„) is related to the probability current density flow from the well to the barrier. 

Using the instanton approach procedure shortly described in the Sections II and III (for the details see [10], [11]), 
in own turn one can find not only eigen values but as well eigen states. The real-valued instanton wave functions arc 
determined by the action W{X\,X) which is counted from the linear turning point X\ 

(X) = A(a)\X - Xi|- 1 / 4 sin( 7 VK(X 1 ,X) + tt/4) , (4.20) 

where the amplitude A(a) of the wave function acquires maximum values at the poles of (4.16) with widths propor- 
tional to T„. We plotted the functions |A(a)| 2 in Fig. 6. 

V. RESONANCE TUNNELING. 

The phenomcnom of the electron resonance tunneling is a familiar one [27] and observed (see e.g. [28], and for 
more recent references also [29]) in semiconducting heterostructures possessing so-called double barrier potentials (see 
Fig. 8). The phenomenom manifests itself as peaks in the tunneling current at voltages near the quasistationary 
states of the potential well. The physical mechanism of the resonance tunneling can be understood as a constructive 
interference between the wave reflected from the left barrier and outgoing to the left from the well. 

In the instanton method the total transmission coefficient T is determined by the second order turning points of 
the double barrier potential (i.e. the maxima of the potential) and according to the procedure shortly described in 
the Section II, T reads as 

| T |2 _ ^lT R (51) 

(1 - V(l + 7rr L )(l + Trr^)) 2 + 4^(1 + 7rr L )(l + ttTr) cos 2 ( 7 W^) ' 



where we designate 



as above (2.12) 



Tl.k = - exp(-7ra L ._ R ) , (5.2) 

IT 



a L . R = 2— i- — '- , (5.3) 

and the action in the classically admissible region analogously to (3.19) is 

7 Vt^ = 7 IU* - <f>{a L ) - cj>{a R ) . (5.4) 
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In the resonance region, where according to the stationary quantization rule 



the transmission coefficient from (5.1) is 



and far from the resonance 



4r r ,r 



Ll R 



4 cos 2 {-)Wr) ' 



(5.5) 



(5.6) 



Thus we found the resonance amplification of the transmission. For the symmetric case in the resonance T = 1, and 
the interference suppresses completely reflection. In the opposite case of strongly asymmetric barriers T from (5.1) 
is almost coincided with the transmission coefficient for the highest barrier, and an influence of the lower barrier is 
suppressed by the interference. In Fig. 9 we show the energy dependence of T for the symmetrical structure of the 
barriers. The resonances become broader when the energy approaches to the potential barriers top, and disappear 
at higher energies (above the top). It is worthwhile to stress that the instanton solution of the resonance tunneling 
problem allows us to study the phenomenon! in a very broad energetical region, including the states near the barriers 
tops. 

Finally let us present the connection matrices for the found above instanton solutions. The corresponding matrix 
can be found as the product of two connection matrices connecting instanton solutions near the second order turning 
points (see (A9) and (3.20)), and the diagonal shift matrix (A6): 



' 2tt exp(7r(a L + a R /A) exp(i7W*){r[(l + ia L )/2]T[{l + ia R )/2}}- 1 + cxp(n(a L + a R )/2)(exp(-i 1 W*) 
-z\/27exp[(37r(c^ + a L )/8] (exp(*7W*)r- 1 [(1 + ia R )/2] + exp^^IF*^" 1 [(1 - ia L )/2]) 



(5.7) 



i\/2^exp[(37r(a fl + a L )/8] (exp^FF^r^Kl + ia L )/2] + expi-i-fW^T' 1 [(1 - ia R )/2}) \ 
27rcxp(7r(a L + a R /A) exp(-i-yW*){T[(l + ia L )/2]T[(l + ia R )/2}}- 1 + cxp(7r(a L + a R )/2){exp(i-fW*) J 

Here as always we designated by W* the action between the turning points (in this case between the second order 
turning points). 



VI. CONCLUSION 



Our paper could be considered as a formal one, in the sense that we asked theoretical questions that most of 
solid state or chemical physics experimentalists would not think to ask. However, the answering of these very basic 
questions can be illuminating. 

Let us sum up the results of our paper. Within the framework of the instanton approach we derived accurate 
analytical solutions for a number of Id scmiclassical problems, and checked numerically the results. As an illustration 
of the method we consider a simple quantum mechanical problem - penetration of a particle through the parabolic 
potential barrier. In this case the instanton solutions which are the asymptotics of the Weber equation are exact. 
The second investigated problem is related to the description of highly excited states in a double - well potential. For 
sake of concreteness and simplicity we study a quartic anharmonic A 4 potential. The instanton approach enables 
us to reproduce accurately not only asymptotic behavior but also a crossover region from single well to double wells 
quantization (in a contrast with the standard WKB approach which fails to describe the crossover behavior). The 
analogous problem for A 3 anharmonic potential is also studied, and the instanton method allows to study resonance 
broadening and collapse phenomena. Besides we investigated so-called resonance tunneling phenomena, interesting 
not only in its own right but as well playing a relevant role in many semiconducting double barrier structures. We 
computed as well the connection matrices, which provide a very efficient method of finding semiclassical solutions to 
the Schrodinger equation in potentials having several turning points. It is also useful for developing good analytical 
approximation. 

All examples selected in our paper to illustrate the efficency of the instanton approach, belong to the fundamental 
problems of chemical dynamics and molecular spectroscopy (see e.g. the monography [9] and references herein). 
Symmetric or slightly asymmetric double-well potentials are characteristic for molecules and Van der Waals complexes 
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with more than one stable configurations [14], [15], [16], [17]. The states of such systems closed to the barrier top 
(theoretically described by the instanton approach in our paper) are not easy to investigate experimentally since 
optical transitions between these states and localized ones, are typically inactive. However just these states are most 
relevant for radiationless evolution of highly excited states. In a certain sense these states have a double nature 
(localized - dclocalized) and the nature manifests itself in the form of wave functions which contain simultaneously 
the both components: localized in one from the wells, and delocalized between the both of the wells. As a consequence 
of it, any initially prepared localized state will evolve via formation and decay of these states. Our calculations are 
intended to pave a way to investigate this class of problems using the wave functions computed within the instanton 
approach. 

The states close to the barrier top of decay potentials govern of thermally activated over-barrier transition ampli- 
tudes. For the low energy states the main reduction factor is the tunneling exponent, while the contribution of the 
highly excited states is limited by the Boltzmann factor. The energetical width of the region dominating in the total 
transition rate is postulated traditionally in the transition rate theory [31] as being of the order of the temperature 
T. However our results presented in the Section III predict another estimation. Instanton calculations demonstrate 
that the intermediate region between the quasistationary r <C £1 and dclocalized states could be much larger than 
T, namely of the order of 0. It means that there is no sharp boundary between quasistationary and delocalized 
states and all of the states within the interval V* — £1 , V* + T equally contribute to the total rate constant for the 
penetration through the barrier. 

One more point should be emphasised. Recently in [18] has been shown that quantum irreversibility phenomena 
occur when the spacing between neighbouring levels of a deeper well becomes smaller than the typical transition 
matrix element. Obviously this criterion can be also applied to the states near the barrier top. Note that for the 
low energy states the asymmetry providing irreversible behavior should be very large, whereas for the states near the 
barrier top, the condition to have the ergodic behavior is not so severe, it is sufficient that the asymmetry of the 
potential is comparable to the barrier height. 

Applications of the method and of the results may concern also the various systems in physics, chemistry and 
biology exhibiting double level behavior and resonance tunneling. Literally speaking in this paper we dealt with the 
microscopic Hamiltonians. However, thanks to the rapid development of electronics and cryogenic technologies, it has 
become possible to apply the same Hamiltonians to study cases where the corresponding variables are macroscopic 
(e.g. magnetic flux through a SQUID ring, or charge or spin density waves phase in a certain one dimensional solids). 
For example in the paper we studied tunneling processes in which a system penetrates into a classically forbidden 
region (a potential barrier). It is an intrinsically quantum effect with no classical counterpart, but nevertheless it can 
take place for macroscopic systems, and the tunneling of a macroscopic variable of a macroscopic system (e.g. spin 
or charge tunneling in atomic condensates trapped in a double- well potential [32]) can be also investigated by our 
method. 

With this background in mind our results are also intended to clarify different subtle aspects of tunneling phenomena. 
An example was given at the end of the Section III , where we found the geometrical phase acquired by a particle 
tunneling through a potential barrier. This phase can be tuned by the particle energy and by the barrier shape, and 
specific interference phenomena might occur. The observation of oscillations related to this geometrical phase in real 
systems has proved challenging. Of course since the forms of the model potentials we used are rather specialized (and 
besides only Id), we cannot discuss the behavior for general cases with full confidence. Nevertheless, we believe that 
the instanton approach employed in the present work should be useful in deriving valuable results for the general and 
multidimensional potentials as well. 

Essentially that in the instanton method, we discussed in the paper, just observing a few classical trajectories suffices 
to develop a qualitative insight for quantum behavior. Though (as we illustrated in a number of particular examples 
considered in the paper) relied upon in this way, semiclassical instanton approach is much more than a qualitative 
picture. As an approximation, the instanton method can be surprisingly precise. Note also that results presented 
here not only interested in their own right (at least in our opinion) but they might be directly tested experimentally 
since there are many systems where the model investigated in the paper is a reasonable model for the reality. 

The theory presented in our paper could be extended in several directions. One very interesting question is how 
our quantum mechanic instanton formulas (e.g. (4.12) - (4.14) for the tunneling rate in the anharmonic X 3 decay 
potential) are modified by interactions with surrounding media (see e.g. [30], where WKB approach was used to study 
the time evolution of quantum tunneling in a thermally fluctuating medium) . However the theoretical modelling of 
this case is hampered by lack of detailed knowledge of the medium and of the interaction with it. A more specific study 
might become appropriate should suitable experimental results become available. A simple criterion for the strength 
of interaction with an environment (or by other words for the effective temperature) where roughly the crossover 
from thermally activated classical to quantum mechanical decay can be found easily by equating the corresponding 
Arrhenius factor and characteristic frequency "oscillations" inside the barrier £1* (see (A10)). 

All of the potentials investigated in our paper can be considered in a number of realistic cases as effectively resulting 
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from avoiding of adiabatic level crossing in the situation, when the adiabatic splitting is so large that any influence of 
the upper adiabatic states on the lower states can be neglected. Certainly in a general case of an arbitrary coupling 
strength, this interaction of higher and lower adiabatic states must be taken into account, and the tunneling matrix 
elements are accompanied by corresponding Franck - Condon factors arising due to violation of Born Oppcnhcimcr 
approximation. We defer these problems to the further, although no any doubt the instanton approach is useful for 
this kind of problems as well. 

Acknowledgments 

The research described in this publication was made possible in part by RFBR Grants 00-03-32938 and 00-02-11785. 

APPENDIX A 

For the standard basic WKB solutions following [2] we will introduce shorthand notations 

(o,z)^( g (z))- 1 / 4 exp( i7 VF(z)), (Al) 

and 

(z, o) = (<z(z))- 1/4 exp(-i 7 W(s)) , (A2) 

The position of the turning point is designated by o and not essential if we are looking for the solutions in the region 
\z\ 3> 1. According to the definitions (2.18), (2.19) on the Stokes lines one should add the dominant solution times a 
certain constant (Stokes constant) to the subdominant (decaying) solution, while on the anti-Stokcs lines the dominant 
and the subdominant solutions are exchanged. To find the Stokes constant we must match the both solutions going 
around the turning point and taking into account the cut on the complex z plane (see Fig. 10). 
Let us consider first the linear turning point 

q(z) = -z, (A3) 

when the classically admissible region corresponds to X > 0. For this case we have 3 Stokes lines, 3 anti-Stokes lines, 
one cut, and therefore 7 different regions on the complex z-plane where the functions (Al), (A2) should be matched, 
and after that 3 Stokes constant should be determined. After not very sophisticated but rather tedious algebra we 
are ending with all three Stokes constants 

Ti = T 2 = T 3 = i , 

and the connection matrix related the coefficients of the combinations of the basic solutions (Al), (A2) in the classically 
are ending (Ai,A 2 ) and in classically admissible {A2,B 2 ) regions are 

where 

M- = exp (-CL\ ( 6X P^/4 { )(%° M 7 i7r /£ ) (A5) 
1 V 4/ \ cxp(-i7r/4) (1/2) exp(i7r/4) J v ; 

Analogously for the other linear turning point q(z) = +z, the connection matrix M + is Hermitian conjugated to M~. 
The variation of the coefficients in the region between two independent linear turning points z\ and z 2 are determined 
by the diagonal matrix L 

L = ( ° XP( 0^ ) eM+ijW*) ) (A6) 

where 



f Z2 

W* = ^qjz)dz. 

J Z\ 
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And finally for the solutions in the classically forbidden regions X < X\ and X > X^ the connection matrix is the 
direct matrix product of the found above matrices 

M = M+LM- . 

To generalize the procedure for the second order turning points we should find the connection matrices related the 
basic solutions to the Weber equation, namely 

(o,z)^(z)"cxp(-z 2 /4), (A7) 

and 

(z,o) = (z)- 1/ - 1 exp(z 2 /4). (A8) 

In this case we have 4 Stokes lines, 4 anti-Stokes lines, one cut, and therefore 9 different regions, where the solutions 
should be matched (see Fig. 11a as an illustration). Four Stokes constants are 

fx = T' 2 _1 (exp(i27ri/) - 1) , T 3 = f* , T 4 = - f 2 cxp(-z2ttz/) . 

And from the known asymptotic of the parabolic cylinder functions one can get the last remaining Stokes constant 

Thus the connection matrix for an isolated second order turning point can be represented as 



— T2 cos(7w) 
cos(7w) — sin 2 (7ri/)/T 2 



(A9) 



For example this depending on the energy e matrix determines instanton semiclassical solutions for the harmonic 
oscillator e = v + 1/2. It can be checked by explicit calculations that for the harmonic oscillator the same form (A9) 
of the connection matrix holds as well for WKB approach. The difference could appear only from anharmonic terms 
in the potential. However for low energy states e/7 <C 1 anharmonic corrections are small and up to the second order 
over these corrections terms the both connection matrices (instanton and WKB) coincide. 

The connection matrix describing the variation of the coefficients at the basic solutions (A7), (A8) between two 
second order turning points X® and X\ for the symmetrical double-well potential is (compare with the analogous 
matrix (A6) for two linear turning points) 

/(n!/V2^)(fio7M)^ +1/2 ex P ( 7 ^) \ 

V (V2^/n!)(r»o7/«*) 1/+1/2 cxp(- 7 T^|)y' 1 ' 

where the instanton action 

x 

W* E = [ ? ' y/2{V{X)-{th)dX t (All) 
Jx° 

and fi* is the characteristic "oscillation" frequency in the barrier (i.e. in the classically forbidden region). For the 
asymmetric double-well potential in the region between the second order and the linear turning points the analogous 
to (A10) matrix is 

(n!/V2^) 1 / 2 (^o 7 /^*)- (1/2)(,y+1/2) exp(7^) \ 

(V2^/n!) 1 /2(fi o7 /^)( 1 /2)(-+i/2) cxp( _ 7M /* )y ) ^ lz > 

All given above matrices allow us to find any other connection matrices we need for all particular examples considered 
in the main text of the paper. Any of them can be constructed as a corresponding product of (A4), (A5), (A6), (A9), 
(A10), (A12). It is worth noting at the very end one general property of the connection matrices, namely that for all 
bound states the connection matrix is real-valued, and for continuum spectrum states, off-diagonal elements of the 
connection matrix are complex. 

Analogously for the problem of tunneling through the potential barrier V{X) = —(1/2)X 2 , all Stokes and anti- 
Stokes lines are turned by the angle 7r/4 (see Fig. lib) with respect to the corresponding lines for the parabolic well 
(V(X) = (1/2)AT 2 considered above, see Fig. 11a). The connection matrix for the tunneling through the barrier 

/ (Si -iexp(7ra/2) \ 

^iexp(7ra/2) 5f 1 (cxp(Tra) + 1) J [ > 
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where a = i(2v + 1) and S\ is the Stokes constant on the first quadrant bisectrix (see Fig. lib). To find the Stokes 
constant Si one has to match the sum of the incident and of the reflected waves to the solutions of the Weber equation 
at X — > — oo and to the transmitted wave at X — > oo. It leads 



Sl = r[(i + L)/2] exp( ™ /4) ' 
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Figure Captions. 

Fig. 1 

The phase of the reflected from the parabolic barrier wave: 

(1) - exact quantum and instanton solutions <po', 

(2) - WKB solution c^ KB ; 

dashed line is the difference <po — (^ KB ■ 
Fig. 2 

The Stokes (solid) and anti-Stokes (dashed) lines for the two real-valued turning points X\^ with the surrounding 
contours I and /'. On the contour II the Stokes lines for the Airy equation asymptotically matched to the lines for 
the Weber equation. The cut is depicted by the wavy line. 

Fig. 3 

The same as Fig. 2 but for the case of the two pure imaginary turning points iX\^- 
Fig. 4 

The dimcnsionless tunneling splitting A/Oo for the anharmonic X A potential near the barrier top: 

(1) - exact quantum and instanton calculations; 

(2) - WKB result. 
Fig. 5 

The X 3 anharmonic decay potential. 
Fig. 6 

The amplitude of the wave function localized in the potential shown in Fig. 5 (the dashed line is a non-resonant 
part of the amplitude, and 7 = 101). 
Fig. 7 

The model two barrier potential structure for the resonance tunneling. 
Fig. 8 

The transmission coefficient for the potential shown in Fig. 7 (7 = 54). 
Fig. 9 

The Stokes (solid) and anti-Stokes (dashed) lines in the vicinity of the linear turning point V(X) = —X. The cut 
is depicted by the wavy line, and the Stokes constant are Ti, T 2 , and T 3 . 
Fig. 10 

The Stokes and anti-Stokes lines in the vicinity of the second order turning points (the same notations as in the 
Fig. 9): 

(a) V(X) = (1/2)X 2 ; 

(b) V(X) = -{l/2)X\ 
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